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Abstract 

We present central-to-peripheral nuclear modification factors, -Rcp, for the pr distri- 
butions of Kg, A, A, and negatively charged particles, measured at central rapidity 
in Pb~Pb collisions at top SPS energy. The data cover the 55% most central fraction 
of the inelastic cross section. The Kg and A Rcp{pt) are similar in shape to those 
measured at ^/Snn = 200 GeV at RHIC, though they are larger in absolute value. 
We have compared our Kg i?cp data to a theoretical calculation. The prediction 
overestimates the data at pT ~ 3-4 GeV/c, unless sizeable parton energy loss is 
included in the calculation. 

Key words: nucleus-nucleus collisions, strange particles 
PACS: 12.38.Mh, 25.75.Nq, 25.75.Dw 



1 Introduction 



The quenching of high transverse momentum (py) particles in central heavy- 
ion collisions is one of the main discoveries at BNL-RHIC. The effect is quan- 
tified using the nuclear modification factor: 

R (n\- ^ iPNSJdpTdy 



where (A^coii)c average number of nucleon-nucleon (NN) collisions for 

nucleus-nucleus (A A) collisions in a given centrality class C. The nuclear mod- 
ification factor would be equal to unity if the AA collision were a mere su- 

^ Corresponding authors: G.E. Bruno, via Orabona 4, 1-70126 Bari (Italy), 
giuseppe.bruno@ba.infn.it; A. Dainese, via Marzolo 8, 1-35131 Padova (Italy), 
andrea . daineseSpd . inf n . it 

^ Permanent address: University of Udine, Udine, Italy 



2 



perposition of A^cou independent nucleon-nucleon collisions. In central Au-Au 
collisions at a c.m.s. energy per nucleon-nucleon pair of ^/s^ = 200 GcV, 
the PHENIX and STAR experiments have measured a suppression by a fac- 
tor 4-5 with respect to unity in it^AA for Pt~5 GeV/c, independent of the 
particle species [1]. A similar suppression at high pr is observed also in the 
central-to-peripheral nuclear modification factor 

p . ^ (^coll)p ^ d'N^JdpTdy 
(A^coii)c d^NlJdprdy 



where a class P of peripheral nucleus-nucleus collisions replaces the pp ref- 
erence (see e.g. [2]). The measured suppression is interpreted as being due 

to energy loss of the hard partons traversing the high-density QCD medium 
expected to be formed in high-energy heavy-ion collisions [3]. Parton energy 
loss would predominantly occur via the mechanism of medium-induced gluon 
radiation [4,5]. 

Preliminary results [6] from the RHIC run at y^s^ = 62.4 GeV show a Raa 
suppression by about a factor 3, not much smaller than that observed at 
^Snn = 200 GeV. This motivates the search for parton energy loss effects at 
even lower c.m.s. energy, i.e. in Pb-Pb collisions with ^s^n = 17.3 GeV at 
CERN-SPS. Possible indications of such effects were found by the WA98 Col- 
laboration [7], who observed a suppression of the tt^ Rcp- In addition, a recent 
71"" Raa analysis [8], using the WA98 Pb-Pb data [7] and a parameterization 
of a wide compilation of pp data at similar energy available in the literature, 
favours a scenario of significant energy loss. 

We investigate the presence of energy loss effects in Pb-Pb collisions at ^ s^n — 
17.3 GeV by measuring Rcp{pt) for Kg, A and A particles, and for unidenti- 
fied negatively charged particles, h~. The relative behaviour of the Rcp{pt) 
patterns for Kg and A particles is also expected to be sensitive to parton coa- 
lescence effects [9] in the hadronization dynamics. It has been suggested that 
these effects occur at RHIC energy (see e.g. [10]). 

In section 2 the experimental setup, the data sets and the analysis procedures 
are described. The 7?cp results are presented in section 3 and they are com- 
pared in section 4 to other experimental results at SPS and RHIC energies as 
well as with theoretical calculations. 



2 Apparatus, data sets and analysis 



The NA57 apparatus, described in detail in [11], was designed to study the 
production of strange and multi-strange particles in fixed-target heavy-ion col- 
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Table 1 

Average number of participants and of NN collisions with their systematic errors, 
as a function of centrality. 



Class ("J 


^0 a.PVPh 

" incl / 


X"' ' part; / 


(A^oii) 

' coil / 


0- 


5.0% 


345.3 ± 1.7 


779.2 ± 26.6 


10.0 


-20.0% 


214.7 ±5.8 


421.7 ±26.1 


20.0- 


-30.0% 


143.0 ± 6.6 


247.7 ±21.5 


30.0- 


-40.0% 


92.6 ±6.4 


140.5 ± 16.2 


40.0- 


-55.0% 


49.5 ± 5.0 


63.8 ± 9.8 



lisions by reconstructing their weak decays into final states containing charged 
particles only. Tracks are reconstructed in the 30 cm-loug silicon telescope: an 
array of pixel detector planes with a cross section of 5 x 5 cm^ placed inside 
a 1.4 T magnetic field normal to the beam direction. The telescope is inclined 
in the non-bending plane by a 40 mrad angle relative to the beam line and it 
points to the target located 60 cm from the first detector plane. The accep- 
tance covers about half a unit in rapidity at central rapidity and transverse 
momentum larger than about 0.5 GeV/c. For the Pb-Pb runs, the centrality 
trigger, based on charged multiplicity, was set so as to select approximately 
the most central 60% of the inelastic collisions. 

The results presented in this letter are based on the analysis of data collected 
during the 1998 (about 2.1 • 10^ events) and 2000 (about 2.5 • 10^ events) 
Pb-Pb runs with beam momentum of 158 A GeV/c. 

The collision centrality is determined using the charged particle multiplicity 
A^ch in the pseudorapidity range 2 < < 4, sampled by the microstrip silicon 
detectors (MSD) as described in [12,13]. Ach is related to the centrality assum- 
ing Ach = q ■ A^part modified Wounded Nucleon model) [12], where Apart is 
the number of participants, i.e. nuclcons participating in the primary nuclcon- 
nucleon collisions, estimated from the Glauber model [14]. The A^h distribu- 
tion is well described by a Wounded Nucleon model fit with a — 1.0 [13]. 
The inelastic cross section extracted from the fit [12], using an inelastic non- 
diffractive NN cross section (x^^^i = 30 mb [15,16], is af^^^^^^ = 7.26 b. Five 
centrality classes arc defined, with Ach hmits corresponding to given fractions 
of CTj^ef^^- eauch class the average values of Apart and AcoU are calculated 
from the Glauber model, with Pb nucleus Woods-Saxon density-profile param- 
eters as given in [17] . In Table 1 we present the definitions of the five centrality 
classes along with the corresponding values of (Apart) and (Acou) with their 
systematic errors, estimated by varying the a parameter in the range [1.0,1.1], 
the Woods-Saxon parameters within their tabulated uncertainties and the in- 
elastic non-diffractive NN cross section at ^/s ~ 17 GeV within its systematic 



4 



'lo 20000 — < Pt < 0-6 GeV/c 

"q _ -J/ 

c\i 1 5000 — 



E 10000 
5000 



2500 

0- 



1.1 

M[p7i"] (GeV/C 



1.4<PT<1.6GeV/c 

> I 
■i 8000- 



« 6000 — 



0) 4000 — 
2000 — 



"•^ L- 



0.4 0.5 0.6 

M[iiV] (GeV/c^) 

1 0000 hO.5 < Pt < 0.6 GeV/c 
CM 7500 — 



^ 5000 — 



1.2 

2i 



8000 — 1 .4 < Pt < 1 .6 GeV/c 

^ : \^ 

<N 6000 — 



H 4000 — 
2000 — 
— 



1.1 

M[p7c"] (GeV/c^; 



> 

(15 

CM 

W 
0) 



0.5 0.6 
M[jiV] (GeV/c^) 



1.2 



- 2.5 < Pt < 3 GeV/c 



200 



100 



i 



| L..tV...,l'.^.l»'" ^ I 

1.4 0.5 



0.6 

2-. 



M[7cV] (GeV/c^) 



> 
<1> 

CM 400 — 
w 

c 
<u 

200 — 



2.5 < Pt < 3 GeV/c 



1.1 



1.2 



M[p7c"] (GeV/c^) 



Fig. 1. Invariant mass distributions for Kg (upper row) and A (lower row) candidates 
in different px ranges. Arrows indicate the nominal masses of the two particles and 
shaded areas indicate the ranges considered for the analysis. No centrality selection 
is applied. 

error estimated to be 1.5% [15,16] ^ . Note, however, that the contribution due 
to the uncertainty on o";^^ is negligible for (iVpart), while for {Nco\\) it ranges 
from 1.5% in the most central class to 1.1% in the most peripheral, but it 
cancels out in the ratio (Ai'coii)p/ (iVcoi^c used for -Rcp- 



Strange particles are reconstructed using their decay channels into charged 
particles: Kg n^ir', A vr~p, and A vr+p. The selection procedure is 
described in detail in [18,19]. The main criteria are the following: (a) the two 
decay tracks are compatible with the hypothesis of having a common origin 
point; (b) the reconstructed secondary vertex is well separated in space from 
the target; (c) the reconstructed candidate points back to the primary vertex 
position. For Kg, the two decay tracks are required to miss the interaction 
vertex by applying a pr-dependent cut on the product of their impact param- 
eters (distances of closest approach to the interaction vertex in the bending 
plane). Ambiguities among Kg, A and A are eliminated by means of cuts [18] 
on the Podolanski-Armenteros plot [20]. The final invariant mass distributions 
for Kg and A candidates in different px ranges are shown in Fig. 1. The shaded 
areas correspond to the windows used for this analysis; the size of these win- 
dows is smaller at low pr, where the invariant mass resolution is better. The 

^ The (A'part) values given in Table 1 for the most central and the most peripheral 
classes differ slightly from those published in [13] due to a refined treatment of the 
trigger effect in the fit of the multiplicity distribution done for the present analysis. 
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Fig. 2. (colour online) Ratio of the average correction weights for a sample of central 
events (0-10%) and a sample of peripheral events (30-55%), as a function of pT and 
of the particle species. 

statistics amounts to 1.8 ■ 10^ Kg, 0.7 ■ 10^ A, and 0.1 ■ 10^ A. We have esti- 
mated and subtracted the combinatorial background using the event-mixing 
technique on a subsample of events representative of the full statistics [19]. 
The background fraction in the selected mass window is negligible for our 
most peripheral class (40-55%) and it increases going to more central classes, 
where it also increases from low to high transverse momentum. In the worst 
case, in class 0-5%, the subtracted background fraction amounts to: (6 ±2)% 
for Kg with pr > 2.5 GeV/c; (10 ± 3)% for A with pr > 2.5 GeV/c; (8 ± 3)% 
for A with pr > 2 GeV/c. 

Negatively charged particles, h~ , are selected from a sample of good-quality 
tracks (clusters in more than 80% of the telescope planes, less than 30% of 
the clusters shared with other tracks) using an impact parameter cut to en- 
sure they come from the interaction vertex. The residual contamination of 
secondary tracks (decay products of weakly-decaying particles) has been es- 
timated to be of about 2%, independent of pt, for the most central class of 
events, and lower for the other classes. The statistics amounts to 10® particles. 

We calculate Rcp{pt), Eq. (2), using px distributions which are unweighted, 
i.e. not corrected for geometrical acceptance and reconstruction/selection ef- 
ficiency. For other analyses, in particular for the measurement of the pt- 
integrated strange particle production yields [21], we adopted a procedure, 
specifically developed for rarer signals like S and Q, in which every selected 
particle was assigned a correction weight, calculated on the basis of a Monte 
Carlo simulation [19]. This correction is very time consuming and, for the 
more abundant signals {h~ , Kg and A), it was calculated only for a repre- 
sentative subsample of the available statistics. For the present study, we have 
verified that these weights do not depend on the event centrality over the full 
transverse momentum range covered (see Fig. 2). Using unweighted spectra 
results in a negligible systematic error on -Rcp(Pr) as compared to the other 
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contributions discussed in the following. 



3 Results 

We use class 40-55% as the reference peripheral class in the denominator of 
it!cp, see Eq. (2), and vary the 'central' class in the numerator from 0-5% 
to 30-40%. We estimated possible residual systematic effects by comparing 
the Rcp results corresponding to the year 1998 sample and to the year 2000 
sample, which were also processed with different versions of the reconstruction 
software. We find compatible results, with point-to-point differences smaller 
than 10% for h~, Kg and A-l-A. We assign 'reconstruction' systematic errors 
of 5% for all four species. For strange particles, a small contribution to the 
systematic errors has been recognized as due to the procedure used to remove 
fake tracks, which yield duplicated candidates, particularly in central collisions 
that have higher track density^ . For K3, A and A, the resulting error, esti- 
mated to be at most 3% at high pt for the ratio 0-5%/40-55% and smaller for 
the other ratios, has been added in quadrature to the 'reconstruction' system- 
atic error and to the systematic error due to background subtraction, which 
is also at most 3%. 

Figure 3 shows our Rcp{pt) results for four different 'central' classes (C). 
The error bars are obtained as a quadratic sum of the statistical errors and 
the pj-dependent systematic errors (< 7%). The shaded bands centered at 
i?CP = 1 represent the pT^-independent systematic error due to the uncertainty 
in the ratio (iVcoii)p/ {Ncoii)c, while the shaded bands at low px represent 
the Rcp values corresponding to A^part-scaling, with the band indicating the 
systematic error due to the uncertainty in the ratio (A^part)c / (-^part)p- We first 
focus on the 0-5%/40-55% ratio. At low pr (~ 0.5 GeV/c), the h', Kg and 
A patterns are compatible with A^part-scaling, while the A points are clearly 
below. As Pt increases, Rcp for the Kg approaches one; Rcp for the A behaves 
differently from the Kg above pr ~ 1 GeV/c, reaching a value of about 1.5, 
as does the A. For pt values below 2 GeV/c, the it!cp for negative particles 
is dominated by negative pions, and stays below the corresponding values for 
Kg. At higher p^ the Rqp lies between those for Kg and A; in this region 
there may be significant K^ and p contributions in the sample of negative 
particles (note that p/tt" ~ 0.8 for Pt ~ 3 GeV/c in central Au-Au coUisions 
at ^Snn = 200 GeV [22]). For less central coUisions, the Kg Rqp exhibit a 
small enhancement in the range > 1.2 GeV/c, while within errors the Rcp 
for the other particles do not vary with respect to 0-5%/40-55% Rcp- 



^ For the h analysis, fake tracks were avoided by imposing tighter track-quality 
conditions. 
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Fig. 3. (colour online) Centrality dependence of Rcp{pt) for h~ , Kg, A and A in 
Pb-Pb collisions at y/s^ =17.3 GeV. Shaded bands centered at Rcp = 1 represent 
the systematic error due to the uncertainty in the ratio of the values of (A'^coii) in 
each class; shaded bands at low pT represent the values expected for scaling with 
the number of participants, together with their systematic error. 



The difference between the A and the A at low pt, already observed in the 
centrality dependence of the pt- integrated yields in Pb-Pb collisions [23], is 
pronounced. On the basis of the symmetry of the apparatus and of the signal 
extraction procedure, we are confident that the effect is not caused by an 
experimental bias. It may be due to a centrality-dependent absorption effect 
of A in a nucleon-rich environment. In this respect, note that at SPS energy 
the A/A ratio is measured to be lower by a factor about 2 in lead-lead with 
respect to sulfur-induced collisions [24] and the A yield per participant is 
measured to be lower by a factor 0.71 ±0.05 in p-Pb [18] with respect to p-Be 
collisions [23]. 
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4 Comparisons and discussion 



In Fig. 4 we compare our results to Rcp measurements at the SPS and at 
RHIC. In the left-hand panel, the WA98 7r° data [7] for the ratio l-6%/22- 
43% in Pb-Pb collisions at ^/s^ = 17.3 GeV are plotted together with the 
NA57 h~ and Kg data for the same centrality classes^ . Using these classes, the 
Kg Rcp is approximately constant at 0.9 for pt > I GeV/c and is significantly 
larger than that measured by the WA98 Collaboration for 7r° (-Rcp ~ 0.6), 
even when taking into account the normalization systematic errors, indepen- 
dent for the two experiments. The h~ data from NA57 are compatible, within 
the systematic errors, with the 7r° data from WA98 for ^ 1.5 GeV/c, where 
the h~ sample is expected to be dominated by vr^. For higher p^, h~ have 
a larger i?cp than vr"; this may be due to increasing contributions from K~ 
and p in the h~ sample. At top RHIC energy, ^/s^ = 200 GeV, the kaon 
Rqp (as measured by PHENIX [25] and STAR [2]) is larger than that of neu- 
tral pions (PHENIX [25]) for Pt ~ 2 GeV/c, while they are similar for higher 
transverse momenta. The observed difference in i?cp between kaons and pions 
at SPS energy, and at RHIC energy for low pt, is reminiscent of the 'Cronin 
enhancement' above A^^cou-scaling originally observed for 3 < < 6 GeV/c 



^ In the WA98 publication [7], the centrality classes are defined using percentiles 
of the measured minimum bias cross section, f^^Agsm b ~ ^ ^- Here, we rename 
the classes in terms of percentiles of the inelastic cross section, c-^gf^^' "^^^ values 
used by the two experiments for the inelastic cross section are very similar: 7.26 b 
for NA57 and 7.41 b for WA98. 
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Fig. 4. (colour online) Left: Rcp{pt) for h and Kg from NA57 and vr^ from WA98 [7] 
in Pb-Pb collisions at ^Snn = 17.3 GeV. Right: Rcp{pt) for Kg and A in Pb-Pb at 
^/s^ = 17.3 GeV (NA57) and in Au-Au at = 200 GeV (STAR) [2]; slightly 

different peripheral classes are employed for the comparison. The bars centered 
at -Rcp = 1 represent the normalization errors; the point-by-point bars are the 
quadratic sum of statistical and systematic errors. 
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in proton-nucleus (pA) collisions at y^s^ values up to 38.8 GeV [26]. This 
enhancement, commonly interpreted as due to initial-state multiple scatter- 
ing (partonic intrinsic transverse momentum broadening), was in fact found 
to increase according to the hierarchy pions-kaons-protons [26]. At RHIC 
energy, owing to the hardening of the pr distributions, the effect would be 
reduced compared to lower energies, as confirmed by preliminary results on 
the particle-species dependence of the d-Au (similar to proton-nucleus) Rcp 
at ^/s^ = 200 GeV [27]. 




The comparison for Kg and A at SPS and RHIC (STAR data for Au-Au at 
^Snn = 200 GeV [2]) is presented in the right-hand panel of Fig. 4. In the Pt 
range covered by our data, up to 4 GeV/c, the relative pattern for Kg and A 
is similar at the two energies, while absolute values are higher at SPS than at 
RHIC, where parton energy loss is believed to have a strong effect. Part of the 
difference between the Rqp values at the two energies may be due to different 
nuclear modification of the parton distribution functions (PDFs). In fact, in 
the ^Bjorkcn range relevant for parton production at a given Pt at the SPS, 
e.g. a; ~ 0.3 for ~ 3 GeV/c, nuclear PDFs are expected to be enhanced 
by about 10-20% (anti-shadowing), while almost no effect is expected for the 
smaller values, x ~ 0.03, relevant for the same px values at RHIC energy 
(see e.g. [28]). At RHIC, the larger Rcp for A with respect to kaons [2] or, 
more generally, for baryons with respect to mesons [25], in the intermediate 
Pt range, 2-4 GeV/c, has been interpreted as due to parton coalescence [10,9] 
in a high-density medium with partonic degrees of freedom. Our data show 
that a similar A-K pattern is present also at ^/s^ =17.3 GeV. We note that 
such a pattern may also be explained in terms of larger Cronin effect for A 
with respect to kaons. 

In Fig. 5 we compare our Kg data to predictions provided by X.N. Wang [29], 
obtained from a perturbative-QCD-based calculation [4,30], including (thick 
line) or excluding (thin line) in-medium parton energy loss. The initial gluonic 
rapidity density of the medium, dNg/dy, was scaled down from that needed 
to describe RHIC data, according to the decrease by about a factor 2 in the 
charged multiplicity per unit of rapidity from RHIC to SPS energy. For the 
0-5%/40-55% -Rcp, the curve without energy loss shows a large enhancement, 
increasing with pt- In the calculation this enhancement arises principally from 
initial-state partonic intrinsic transverse momentum broadening, which is as- 
sumed to be proportional to the number of scatterings that the two colliding 
partons suffer inside the nuclei before the hard scattering, and thus larger for 
central than for peripheral collisions [30]. The magnitude of the broadening is 
tuned on the basis of the original Cronin effect data [26]. A 'Cronin-like' en- 
hancement is observed in the d-Au central-to-peripheral nuclear modification 
factor for charged kaons at RHIC energy [27], while it is not present in our 
nucleus-nucleus Kg data. The curve that includes energy loss, scaled down 
from RHIC as explained above, describes the data better. Moving to less cen- 
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Fig. 5. (colour online) -Rcp(pt) for Kg in Pb-Pb collisions at ^/s^ = 17.3 GeV, 
compared to predictions [29,31] with and without the effect of parton energy loss 
(details in the text). 

tral collisions the predictions with and without energy loss get closer to one 
another, and both are compatible with the data. As a cross-check, we compare 
the i?cp value predicted by X.N. Wang with energy loss to the prediction of 
an independent model of parton energy loss, the Parton Quenching Model 
(PQM), based on the BDMPS formalism [5], that describes several energy- 
loss- related observables at RHIC energies [31]. Since the PQM model does not 
include initial-state effects, the predicted Rcp was rescaled using the Wang 
baseline without energy loss (thin line). The PQM result has an uncertainty 
due to the fact that the medium-induced energy loss becomes of the order of 
the initial parton energy [31]; since this uncertainty is larger for low-energy 
partons, the calculation result is meaningful for pt ~ 4 GeV/c only. The PQM 
result, shown in Fig. 5 by the hatched areas, agrees with the Wang calculation 
with energy loss (thick line). The two models predict a similar energy loss effect 
at SPS energy, i.e. a reduction of the 0-5%/40-55% Rqp (for pt~4 GeV/c) 
by about a factor 2, with respect to the value calculated without energy loss. 



5 Conclusions 



Central-to-peripheral nuclear modification factors for Kg, A, A and in 
Pb-Pb collisions at top SPS energy have been measured as a function of 
Pt up to about 4 GeV/c. At low px, Rcp agrees with A'part scaling for all the 
particles under consideration, except the A, for which the yields at low pr are 
found to increase slower than the number of participants. For 1 GeV/c, 
Kg, A and A show a pattern similar to that observed in Au-Au collisions at 
top RHIC energy, although the Rcp values are found to be larger at SPS. 
At RHIC, this pattern has been interpreted in the framework of models that 
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combine parton energy loss with hadronization via coalescence, at intermedi- 
ate pt, and via fragmentation, at higher p^. The measured Kg 0-5%/40-55% 
Rqp is not reproduced by a theoretical calculation that includes only initial- 
state nuclear effects. The data can be better described by including final-state 
parton energy loss as predicted for SPS energy on the basis of RHIC data. 
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